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Abstract. We discuss the proper definition of the spin current operator in Heisenberg magnets subject 
to inhomogeneous magnetic fields. We argue that only the component of the naive "current operator" 
JijSi x Sj in the plane spanned by the local order parameters (Si) and (Sj) is related to real transport 
of magnetization. Within a mean field approximation or in the classical ground state the spin current 
therefore vanishes. Thus, finite spin currents are a direct manifestation of quantum correlations in the 
system. 



PACS. 75.10.Jm Quantized spin models 
Persistent currents 



1 Introduction 



In a recent Letter pQ and a subsequent paper [2] we 
have calculated the persistent spin currents in meso- 
scopic Heisenberg rings subject to inhomogeneous mag- 
netic fields. We have emphasized the close analogy be- 
tween this phenomenon and the well known persistent 
charge currents in mesoscopic metal rings pierced by an 
Aharonov-Bohm flux. In the ensuing discussions with sev- 
eral colleagues we have become aware of the fact that the 
definition of the spin current operator in Heisenberg mag- 
nets subject to inhomogeneous magnetic fields is not ob- 
vious. In this note we shall attempt to clarify this point. 

A related problem, which will not be discussed in this 
work, is the definition of the spin current operator in semi- 
conducting electronic systems with strong spin-orbit inter- 
actions. Recently, Rashba [3] pointed out that also for this 
case the precise meaning of the concept of a spin current 
is rather subtle. In particular, he emphasizes that spin 
currents in thermodynamic equilibrium, which arise with 
the standard definition of the spin-current operator used 
in the literature, are unphysical and should be regarded 
as background currents which do not correspond to real 
transport of magnetization. A clear understanding of this 
concept is essential for the highly active field of informa- 
tion processing using spin degrees of freedom subsumed 
under the name of spintronics 0|. 

For itinerant systems the spin is an intrinsic property 
of the charge carriers and is carried around with their mo- 
tion. For localized spin systems considered here, transport 
of spin is a consequence of the time evolution of the mag- 
netization. For special cases the transport can be ascribed 



75.10.Pq Spin chain models - 75.30.Ds Spin waves - 73.23.Ra 



to the movement of quasi-particles as magnons or spinons 
and again a simple physical picture emerges 

In this context, it is also interesting to note that in ef- 
fective low-energy models for ferromagnets involving only 
the spin degrees of freedom even the concept of the lin- 
ear momentum is not well defined [B]. In general, the dy- 
namical equation for the spin degrees of freedom have to 
be supplemented by kinetic equations for the underlying 
fermionic excitations. 

At the heart of the ambiguities involved in defining 
a spin current operator both for itinerant systems with 
spin orbit interaction as well as for Heisenberg magnets in 
inhomogeneous fields is the fact that the magnetization is 
not strictly conserved for these systems. Still, the intuitive 
concept of magnetization transport should also be useful 
for these systems and one is therefore led to define effective 
current operators, as we will do in this note for the case of 
a Heisenberg magnet in an inhomogeneous magnetic field. 



2 Problems with the naive definition of the 
spin current operator 

Consider a general Heisenberg magnet with Hamiltonian 



H — — ^ JijSi ■ Sj — g[iB ^ B{- Si, 



(1) 



where the sums are over all sites r< of a chain with pe- 
riodic boundary conditions, the J,j are general exchange 
couplings, and Si are spin-S" operators normalized such 
that Sf = S(S + 1). The last term in Eq. |TTJ is the Zee- 
man energy associated with an inhomogeneous magnetic 
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Fig. 1. Classical spin configuration rhi of a nearest-neighbor 
ferromagnetic Heisenberg ring in a radial magnetic field Bi. 

field Bi = B(ri). We assume that the magnetic field at 
each lattice site is sufficiently strong to induce permanent 
magnetic dipole moments rrii = g^(Si), not necessarily 
parallel to Bi, where (...) denotes the usual thermal av- 
erage. The simplest geometry is a ferromagnetic ring in a 
crown-shaped magnetic field, as illustrated in Fig.fJ This 
geometry is used in the following for illustrative purposes, 
but our arguments are not restricted to this case. The 
Hamiltonian (Tj^) implies the equation of motion 

dS- 

h-^+hiXSi + ^JijSiXS^O, (2) 

3 

where hi = g^Bi. Note that the last term in Eq. (0 can 
be written in the form Y) . Ii^j , with 

li^j — JijSi X Sj . (3) 

By analogy with the discrete lattice version of the equa- 
tion of continuity for charge currents, it is tempting to 
identify Ii^j with the operator whose expectation value 
gives the spin current from site i to site j. In this work 
we shall argue that this identification is only correct for a 
strong homogeneous magnetic field, where in equilibrium 
the expectation values (Si) of the spins at all sites are 
aligned along the same spatially constant direction of the 
field. Then hi x (Si) = 0. Using the fact that equilibrium 
averages are time independent, jjj(Si) = 0, we conclude 
from the equation of motion J5J that the lattice divergence 
of the spin current in the presence of a homogeneous mag- 
netic field vanishes 

£<W=0. (4) 

3 

For a one-dimensional ring with nearest neighbor hopping 
this implies for each site i 

(ii-i+i) + (Ii^i-i) = (Ii^i+i) - (/i-i-i) - , (5) 

so that the same spin current (I) = (_Z^,; + i) flows through 
each link of the ring. However, the equation of motion con- 
tains only the divergence of the current, so that it does 
not fix the value of (I). From the point of view of ele- 
mentary vector analysis this is a consequence of the fact 
that both the divergence and the curl are necessary to 



uniquely specify a vector field. Because the equation of 
motion contains only the divergence, circulating spin cur- 
rents cannot be calculated using the equation of motion. 
In fact, even the definition of the spin current operator in 
a geometry permitting circulating spin currents cannot be 
deduced from the equation of motion. Of course, for a ring 
with a collinear spin configuration we know that (I) = 0, 
so that there are no circulating currents. 

The case of a non-uniform magnetic field is more in- 
teresting. In general, the spin configuration in the ground 
state is then also inhomogeneous. For example, let us con- 
sider a radial magnetic field Bi = \B\ri/\ri\ situated at 
constant latitude i3i — as shown in Fig. ^ We assume 
that the direction rhi = mi/\mi\ of the magnetic mo- 
ments rrii = g/iB(Si) trace out a finite solid angle fl on 
the unit sphere in order-parameter space as we move once 
around the ring. If we consider a nearest neighbor Heisen- 
berg ferromagnet with ~ — J(^ J+ i + 6i_j-i) then for 
\h\ = i?/ib|-B| > JS(2tt/N) 2 the classical ground-state 
configuration rhi is radial as well, with a slightly different 
latitude *& m satisfying pQ 

sin(i? m -■&) = -(JS/\h\) [1 - cos(2tt/A0] sin(2tf m ) . (6) 

The main point of this work is that in the presence 
of an inhomogeneous magnetic field the spin current op- 
erator is not simply given by Eq. ®. The fact that the 
expectation value of Eq. @ cannot be a spin current is 
perhaps most obvious if we consider the simple case of 
classical spins in a star-shaped magnetic field, correspond- 
ing to t? m = 1? = 7r/2 in Fig. ^ In this case Eq. © gives 
for a ring with evenly spaced sites at zero temperature 

Ii^j = J y e z sin(27r/JV) , (7) 

where e z is a unit vector perpendicular to the plane of 
the ring. Note that at the classical level the statics and 
dynamics of a Heisenberg magnet are completely decou- 
pled. Because the classical ground state does not have any 
intrinsic dynamics, it does not make any sense to associate 
a spin current with it which would correspond to moving 
magnetic moments. Furthermore, if the classical Heisen- 
berg model is provided with Poisson bracket dynamics, the 
classical ground state yields a stationary solution, since it 
minimizes the energy. Clearly, such a completely station- 
ary state cannot be used to transport magnetization. We 
conclude that for twisted spin configurations Eq. J2J is 
not a physically meaningful definition of the spin current 
operator. 

3 Effective spin currents with correct classical 
limit 

In order to arrive at a better definition, consider a non- 
equilibrium situation, i.e. start with a given density matrix 
at time t = and let the system evolve according to the 
unitary dynamics generated by the Hamiltonian in Eq. (TJQ) . 
The equation of motion J5J then directly translates to a 



Florian Schiitz et al.: What are spin currents in Heisenberg magnets? 



3 



relation for the local and instantaneous order parameter 
dtiSJt + hi x (Si) t +5^(J <w ) t = 0. (8) 

3 

Here, (...)* denotes an average with respect to the time 
dependent density matrix. It is then reasonable to demand 
that a transport current can lead to an accumulation of 
magnetization, i.e. a change in the magnitude of the local 
order parameter in time. For this magnitude, we obtain 
the equation of motion 

d t \(S i )t\ + J2™i( t )-( I i^3)t = 0, (9) 

3 

where rhi(t) = {Si) t /\(Si) t \ is the time dependent direc- 
tion of the order parameter. Note that only the longi- 
tudinal component of the naive "current operator" ap- 
pears in this continuity equation without source terms. It 
is precisely this contribution which we have identified as 
the dominant one in our spin- wave calculation in £Q . The 
transverse components lead to a change in the direction 
of the local order parameter, but they are largely counter- 
acted by the magnetic field term that acts as a source and 
generates a precession. If one wants to discuss the elec- 
tric fields generated by the magnetization dynamics, one 
cither has to take into account both the current Ii— >j and 
the local precessional motion, or devise a way to make the 
cancelation explicit by including part of the "transverse 
current" in an effective magnetic field. We will attempt 
the second route here. 

That this is a sensible way to proceed can be appre- 
ciated by a simple approximate calculation. In the classi- 
cal ground state, the magnetization aligns parallel to the 
sum of the external and the exchange field. A necessary 
condition for the minimum of the classical energy is the 
invariance under small variations of the directions of the 
magnetization. This leads to the condition fjQ 

hf x {Si) = , hf = hi - J io ( s 3) ■ (10) 

3 

Note that the effective magnetic field contains a part of 
the exchange interaction, which therefore should not be 
included into the definition of the spin current operator 
to avoid double counting. Rewriting the exact equation 
of motion (TJJl in terms of the effective magnetic field hf 
defined in Eq. (|1U|) . we obtain 

3 

where 

If^ = JijSi x [Si - (Sj)] . (12) 

Obviously, 

(i£ 3 .) - Ja [(Si x S^ - (Si) x (Si)} , (13) 

which vanishes identically in the classical ground state, or 
if the spins are treated within the mean-field approxima- 
tion, where the spin correlator is factorized. Physically, 



this is due to the fact that within the mean-field approx- 
imation the Heisenberg exchange interaction is replaced 
by an effective magnetic field, so that the different sites 
are uncorrelated and there are no degrees of freedom to 
transfer magnetization between them. In this work we dis- 
cuss only localized spin models, so that charge degrees are 
not available to transfer magnetization between different 
sites. 



4 New definition of the spin current operator 

The definition of ifj in Eq. (|12[> has the disadvantage 
of not being antisymmetric with respect to the exchange 
of the site labels, although its expectation value is obvi- 
ously antisymmetric. In order to cure this problem and to 
generalize the concept of an effective current operator be- 
yond the mean-field description, we propose the following 
definition, 

Ii— >j = li^j ~ lijilij ' ' (14) 

with the unit vector 



Thus, we interpret only the projection of /«— >j onto the 
plane spanned by the two local order parameters and 
rrtj as a physical transport current. The contribution sub- 
tracted in Eq. I|14l) can be incorporated in an effective 
magnetic field. More precisely, the equilibrium expectation 
value of the exact equation of motion (T2J) can be rewritten 
as 

hfx{S i ) + Y i (i ir + j )=0, (16) 

3 

with the effective magnetic field now defined as 

This reduces to Eq. (|l(Jfl for the classical ground state or 
at the mean-field level, where the correlation function in 
the numerator is factorized. The spin current operator de- 
fined in Eq. i|14|) is manifestly antisymmetric under the 
exchange of the labels, as it should be. It implicitly de- 
pends on the spin configuration via the unit vector 7^, 
so that in twisted spin configurations the spin current op- 
erator is a rather complicated functional of the exchange 
couplings. The fact that the current operator of an inter- 
acting many body system is a complicated functional of 
the interaction is well known from the theory of interacting 
Fermi systems 7 . In particular, when the effective inter- 
action does not involve densities only the construction of 
the current operator is not straightforward [Hj- 

For explicit calculations we use a representation of 
Ii—>i in terms of spin operators quantized in local ref- 
erence frames with the z-axes pointing along rhi, i.e. we 
decompose 

S t = S ? e "> S? = e?-Si, (18) 

q=1,2,3 
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with ef = rhi. One still has a freedom in the orientation 
of the transverse basis {e*,ef}, which can elegantly be 
parametrized, if one uses spherical basis vectors ef = e}± 
ief. We can then write 

e+ = e— ef, (19) 

where {ej,ef} is the special transverse basis where ef = 
= jij. With this notation we obtain the following ex- 
pression for the spin current operator 



_ q- q-„i(ui^j+^j^i) _ 

+S!sre i "^(j ij xrhi) 
-STSje^-'faxrhj)- 



2 



H.c. 



(20) 



where S { 



S} ± iSf = ef ■ St are the usual ladder op- 
erators and sf = Sf = rhi ■ Si. The third and fourth 
terms in this expression couple longitudinal and trans- 
verse degrees of freedom and therefore do not contribute 
to leading order in a spin-wave calculation. The first and 
second summand are dominant for ferromagnetic and anti- 
ferromagnetic rings respectively and have been discussed 
in detail in pQ and 0. For a magnetic field that varies 
smoothly as one moves through the system, the magnetic 
moments on neighboring lattice sites are almost collinear, 
so that in both cases the component of the naive "current 
operator" Jj_>j along the local order parameter is the one 
that really corresponds to the transport of magnetization. 
In pP we had come to the same conclusion by invoking 
the gauge freedom in the choice of the transverse axes of 
quantization. Rotating the coordinate frame around rhi 
corresponds to the gauge transformation 



St -> Sfe ±iai 



(21) 



By this gauge freedom one is then let to identify the 
derivative of the Hamiltonian with respect to the gauge 
field LOi^j as the relevant current operator. A comparison 
with Eq. (121 )[l shows that this is indeed the longitudinal 
component of Ii—>j. A more general gauge invariant for- 
mulation of the Heisenberg model is discussed in |9| (see 
also II) i. In these works, an 0(3) gauge field Ai^.j was 
introduced in a rather formal manner to write the Heisen- 
berg model in a gauge invariant way and to obtain the spin 
stiffness tensor by means of differentiation with respect to 
the gauge field [TT] . 

Note that the procedure adopted in this section is not 
restricted to the isotropic Heisenberg interaction of the 
Hamiltonian (Tj^). For a general bilinear spin-spin interac- 
tion of the form 

— S ' Si ■ JijSj ^ hi ■ Si , (22) 



H 



Jji , the equation 



where Jij is now a 3x3 matrix with Jfj 
of motion remains valid, if the expression for Ii—>j 
replaced by 



With this notation, Eqs. I|14I17|I still hold (provided Jij is 
replaced by the matrix Jy ) and part of the naive current 
operator Jj_ >j can again be absorbed into the definition of 
an effective magnetic field. 



5 Conclusion 

Let us emphasize again that our main point is rather sim- 
ple: The microscopic equation of motion © contains only 
the (lattice) divergence of the spin current operator, which 
is not sufficient to fix its rotational part. A certain part 
of the operator Y]^ Jij Si x Sj leads to a renormalization 
of the effective magnetic field and therefore should not 
be included into the definition of the spin current oper- 
ator, see Eqs. (I10H17II . The physical spin current, which 
corresponds to the motion of magnetic dipoles, must be 
defined such that a purely static twist in the ground state 
spin configuration of a classical Heisenberg magnet is con- 
sidered to be a renormalization of the effective magnetic 
field, and does not contribute to the spin current. To fur- 
ther substantiate our proposal for an effective current op- 
erator, it would certainly be insightful to look for a more 
microscopic derivation by starting from a model involving 
charge degrees of freedom. It would also be instructive to 
explicitly investigate non-equilibrium situations with time 
dependent magnetizations. 

We thank P. Bruno and K. Saito for discussions and for 
urging us to clarify the proper definition of the spin cur- 
rent operator in magnetic insulators. This work was sup- 
ported by the DFG via Forschergruppe FOR 412, Project 
No. KO 1442/5-3. 
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